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Abstract 

We derive the effective angular momentum operator to 1/m^ and one- loop order in non- 
relativistic quantum electrodynamics (NRQED). In both dimensional and three-momentum-cutoff 
regularization schemes, we obtain the non-relativistic expansion for the spin and orbital angular 
momentum operators of the electron and photon, respectively. Our result is useful in precision 
calculations of the angular momentum properties for non-relativistic QED bound states, such as 
atom systems. 
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I. INTRODUCTION 



The spin structure of bound states is a topic of considerable interest in recent years [jj. 
The so-called "proton spin crisis" was generated from the EMC experiment [2| in which it 
was found only a small fraction of the total proton spin is carried by the quark spin. In 
Ref . [sj , a gauge-invariant decomposition of the nucleon angular momentum into the quark 
spin, orbital and the gluon contributions has been introduced. The connection between 
off-forward (or generalized) parton distributions and the quark/gluon angular momenta 
within the proton has been established. The total angular momentum carried by quarks is 
measurable through the deeply- virtual Compton scattering (DVCS) j^j. A large amount of 
works have also been done in studying the angular momentum carried by the gluons, AG, 
and several experiments have been designed to measure this quantity and interesting physics 
has been learnt j^. 

While much of the attention has been focused on relativistic theory, it is illuminating to 
study angular momentum in the non-relativistic limit of QED/QCD. A clear understanding 
of the non-relativistic angular momentum structure is useful for many systems, including the 
hydrogen-like atoms ^ and the heavy quarkonium 0] bound states. The fraction of angular 
momentum carried by the photon in hydrogen atom is relevant to the gluon contribution 
to the spin of the proton. Since the hydrogen atom is a non-relativistic bound state, the 
naive perturbation theory of relativistic QED is not applicable because multiple energy 
scales present. For the ground state, the electron has orbital velocity v ~ where 
ctg^ ~ 1/137 is the electromagnetic fine-structure constant, therefore the loop expansion in 
a^^/v does not converge jsf. As another example, in production of electron-positron pairs 
near the threshold, the interactions are always accompanied with the famous Sommerfeld 
enhancement proportional to a^^^/v. Attempts to calculate these corrections within the 
relativistic theory is difficult analytically, usually involving numerical solution of the Bethe- 
Salpeter equations. 

Effective field theory approach sheds new light on the problems described above. In recent 
years, effective theories including non-relativistic QED (NRQED) (of and the non-relativistic 
QCD (NRQCD), has become a usual tool in solving non-relativistic bound states. NRQED 
has been used to calculate the hyper fine splitting and lamb shift of the hydrogen system 



with considerable simplification |10l-|12|. NRQCD has been used in analyzing the heavy 
quarkonium inclusive production in colliders [l3| and precision bound-state calculations on 
the lattice. By observing that the Coulomb interactions exchanging momentum at much 
lower scale compared to the fermion mass m, one can integrate out the large momentum 
scale, resulting in higher dimensional effective operators suppressed by powers oi p/m <C 
1. In effective theories, relativistic and radiative corrections to bound state problems can 
be systematically expanded in terms of p/m and Og,^. The non-relativistic bound state 
wavefunction properly sums up all order corrections due to the Coulomb photon exchange. 
Therefore, the NRQED as an effective field theory, is capable of describing the hydrogen 
atom up to any definite order with desired precision in a language familiar in quantum 
mechanics. 

In the non-relativistic limit, the orbital angular momentum L and spin S of the electron 
are both conserved and can be used to classify the energy eigenstates. Moreover, the con- 
tribution from the gauge potential in the orbital part is proportional to the velocity and 
hence negligible. L is then the usual non-relativistic angular momentum. When taking into 
account the relativistic effects and quantum corrections, the angular momentum operator 
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becomes more complicated. For example, for relativistic Dirac Hamiltonian, neither L nor 
S maintains as a good quantum due to corrections of order v^/c^. These are relativistic 
corrections starting at order 0(a^ ). Radiative corrections will further contribute at order 

I 1 ^ EM' 

^Km) H- ^ 

In this paper, we extend the discussions in [ij] and present relevant details of deriving 
the effective angular momentum operator in NRQED. For this purpose, we will construct a 
set of gauge-invariant effective operators in the non-relativistic theory. The matrix elements 
of the spin (orbital) angular momentum for the electron and the gauge fields are calculated 
in full QED and matched to the effective theory up to order a^j^/m^. The main results 
of this paper are Eqs. (fT3|) -( IT5|) and Eqs. (1461) . (!54l) . (1561) . Using the power counting in 
NRQED [nl, [isj, we have applied the one-loop matching result to calculate the angular 



momentum carried by radiative photons in the hydrogen atom up to 0{a^^) in Ref. [14 
The further extension to the angular momentum decomposition in QCD/NRQCD will be 
presented elsewhere. 

The NRQED Lagrangian is first given in Ref. [9|. In Ref. [l^, a power counting rule has 
been established. 



'2m Sm? 2m Sm? Sm? 
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here C^^-pcrmi represents all the four-fermi terms in the Lagrangian and will not be used in 



this paper. The Wilson coefficients in the effective Lagrangian are well known [18|. Up to 
C(aEM) order, 

C2 = C4 = 1; CF = l + ^- c,, = l-^ln4; = 1 + 

rfi = l + ^ln4; d, = ^. (2) 

In the calculation, we will use the "old-fashioned" NRQED [l6|, which is equivalent to the 
intermediate effective theory in pNRQED [13]. In performing the matching, 1/m expansion 
of the matrix elements (so are the effective operators) as described in Ref. [18| is understood. 



It is proved in Ref. pj| that this agrees with the multi-pole expansion method in Ref. [20 
More explanations will follow along with the detailed calculations. 



II. ANGULAR MOMENTUM OPERATORS IN QED AND NRQED AND MATCH- 
ING CONDITIONS 

The question we are interested in is how the total angular momentum is distributed among 
various components, e.g., that carried by electron and photon, respectively. The effective 
non-relativistic angular momentum operators are useful in studying many non-relativistic 
bound state systems, including the hydrogen-like atoms and heavy quarkonia. The purpose 
of the following sections will be the construction of gauge-invariant angular momentum 
operators in NRQED. We start by first reviewing the angular momentum operator in full 
QED. We then use the Foldy-Wouthysen (FW) transformation to get the corresponding 
NRQED operators at tree level. We then use the Noether current method to derive the 
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total NRQED angular momentum operator from the effective NRQED lagrangian. This 
method is accurate to all loops. However, it does not give us the results for individual 
relativistic angular momentum components. Therefore, we write down the Wilson expansion 
of these angular momentum components in terms of the non-relativistic angular momentum 
operators. In the following sections, we determine the Wilson coefficients using the matching 
method. 



A. Angular Momentum in Full QED 

It is straightforward to derive the total angular momentum operator from the QED 
Lagrangian using the Noether current method. Starting from 

>Cqed = —F^^'F,, + ^{ilp - m)^ , (3) 



the conserved QED angular momentum is [21 



Jqed = j d^x j^^y^ + ^^(x x tt)^ + x x (E x B)| = S, + L, + , (4) 



in which 



= / rf^x X X (E X B) . (5) 



The individual relativistic operators Sg,Lg, are gauge invariant and are regarded as the 
electron spin, electron orbital angular momentum and photon angular momentum, respec- 
tively. Alternative separation can be achieved in some particular gauge and frame choices. 

Our goal is to construct the non-relativistic counterparts to Sg, L^, and J^. The result 
within NRQED is gauge invariant, with the manifest separation of high- and low-scale 
physics. 



B. Non-Relativistic Reduction: Foldy-Wouthysen Transformation 

In this subsection, we apply the Foldy-Wouthysen (FW) transformation to the QED op- 
erators in Eq. ([5]). We shall keep in mind that the FW transformation does not yield any in- 
formation about radiative corrections, due to its quantum mechanical nature. Consequently, 
one can obtain only the leading-order Wilson coefficients by transforming the relativistic op- 
erators Sq, hq, and J^, respectively. Any Wilson coefficient starting from order 0(0;^^) or 
higher would not show up in the transformation. Furthermore, the FW transformation only 
rotate the electron wavefunction, leaving the photon sector unchanged. 

Here is a quick review of the FW transformation. The Dirac Hamiltonian describes the 
relativistic fermion containing both positive and negative frequencies. At lower energy, the 
negative frequency part (positron) is subdominant. For an electron state, the physics is 
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dominated by the interactions of the electron itself, and the correction from the positron 
is suppressed by the mass. The latter is reflected by the off-block-diagonal elements in the 
Hamiltonian. The FW transformation utilizes a unitary transformation to block-diagonalize 
the Dirac Hamiltonian up to order 0(m~^). The unitary transformation matrix is 



U 



where in the Dirac representation 





(3cy. ■ TT 




"7 ■ TT" 


exp 


2m 


= exp 


. 2m . 



(6) 



f3 



'1 




"Oct" 


-1 


, a = 


cr 



(7) 



One can use U to eliminate the lower components in a four-component Dirac spinor of 
electron up to order 0{m~^) and a quantum mechanics operator O transforms accordingly 







0(m" 



u{p) -uip) = 

The operators in Eq. ([5]) take the following form under the FW transformation 



(8) 



2 



d^x \I/^x X 7r\I/ -> 



cr 

2" 

X X 77 



[(cr X tt) X TT — TT X (cr X tt)] 



eB 

4m^ 



^ [(cr X tt) X TT — TT X (cr X 7r)l — 1 ijj 

2 LV / \ n J 



5m 



+ / d^x^ {x X [B X V(?/^V)] + X X ^/^UB X (<T X V)l ?/^| , 
J 8m^ 

ci^x X X (E X B) ^ j d^x x x (E x B) , (9) 

where ijj is the non-relativistic electron field and ^ is defined as: 

= 4(7r^p) + (ttV^p') = i^l' (P + P' - 2eA) . (10) 

In the above results, the next-to-leading order relativistic corrections has been properly 
taken into account. In Ref. it has been shown that the relativistic correction to the 
orbital angular momentum of ground state hydrogen atom can be calculated with the above 
leading order operators. As we will see in the following section, these results are also 
consistent with those derived from matching between QED and NRQED at (9(aj,^/m^) 
order. 



C. Total Angular Momentum from NRQED Lagrangian 

In this subsection, we will first derive Jnrqed using the Noether current method. The 
NRQED Lagrangian contains a second-order-derivative term ^F^J^F^^ and calls for spe- 
cial treatments. In Appendix A, we present the generalized formalism to derive the equation 
of motion and the conserved currents, including the energy-momentum current T^^ and the 
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angular momentum current Ai^'^'^, from a general Lagrangian C = 0, 9^0, d^dycj)). For 
NRQED, the equation of motion with regard to field is 

- eV^t^ _ ^ V2(^t^,) + !^ V ■ [^t X D) + - = 0.(11) 

and the total angular momentum 

d=^x I^Z-t (^^j + ^/;t (x X tt) ^ + X [B X V(^AV)] 

^ X X [B X (o- X ^)] ^ + di X X (E X B) 



J 



NRQED 



+ 



Ad. 



-X X (DE X B) + E''(x X V)E'' - B'^(x x V)B" + ExE-BxB 



:i2) 



We will not consider four-fermion operators in this paper. The total angular momentum 
operator Jnrqed in NRQED is manifestly gauge invariant. However, it is difficult to see a 
clear correspondence with the individual relativistic operators. Therefore, we shall use the 
matching method to derive the non-relativistic expansion of the relativistic operators. 



D. Relativistic Angular Momentum Operator in NRQED 

We observe that the effective Lagrangian in Eq. ([1]) still possesses the rotational symmetry 
in 3D, therefore implies the existence of conserved total angular momentum Jnrqed 0. In 
this paper, we will derive non-relativistic effective angular momentum operators up to order 
O (Og^/m^), through order-by-order matching to the full QED. 

The general form of angular momentum operators in NRQED can be constructed as 
follows. First, we write down all the axial- vector Hermitian operators respecting gauge and 
3D rotational invariance. Second, we require these operators contain no time derivatives on 
the fermion fields. The removal of can be done by field redefinition as shown in (isl ]. 
Any operator containing can be rewritten as linear combinations of other operators by 
the equation of motion without changing the on-shell matrix elements. Third, there will be 
non-local operators in Lg and J^, containing the space coordinate x explicitly. 

Here we write down the general form of the effective operators, with the Wilson coefficients 
to be determined. For electron spin. 



off 



/ dPx I a^i/j'^—i/j + -^^ifj"^ [(cr X tt) X 77 — TT X (cr X tt)] ip + —^ip^'Bi/j 
J I 2 8m-^ Am^ 





Q71 


E X E - B X B 


+ 











m 



[V-E)B- (B- V)E 



0{m~ 



(13) 



The authors in Ref. [2^ have noticed this as weh, and derived the angular momentum in the NRQED up 
to order O ia^ I . 

V EM ' / 



6 



for electron orbital momentum, 



^ 



cff 



J d^x ^da^^^^J + ^^^^ [(''■ X7r)x7r — 7rx(crx tt)] ^ 
+dRi;^{^ X 7r)i; + ^^^B^ + |%x x [B x 



e d.c 



-X X 



edE 
Am 

d. 
m 



[B 



X (cr X 



e d' 

^ + — |x X t/^Uo- X (B X tp 
8m^ 

d^ 



^72 
i2 



X X (^cr X E) -0 + / (i'^x X X (E X B 



E X E-B X B 



V-E B- B - V E 



m^ 



m^ 



E^fx X V)E" - B'^fx X VIB'^ 




(14) 



and for photon angular momentum operator, 



off 



J7 ^ 



[ d^X \ fcrljj^—lp + -^^^^ [{a X TZ) X TT — TT X {(T X tt)] 1p 

J I 2 8m^ 

+/h^/'Hx X 7r)V; + ^^^BiP + ^x X [B X V(0V)] 



efs 



4m 



e f 

X X ?/;UB X (cr X + ^x x ?/;Ucr x (B x 

8m^ 

X X (o- X E) ^1 + j d^xf^y. x (E x B) + 0{m~^) . 



(15) 



In what follows, we will perform the matching of various matrix elements between QED and 
NRQED, to determine all the Wilson coefficients in Eq. ffT3|) -f lT5|) . 
The matching condition is 



(J(/i)) 



QED 



NRQED 



(16) 



Although the primary application of our NRQED effective operators involves bound states, 
we are not obliged to use a bound state to do the matching. For local operators which do 
not contain space coordinate x explicitly, we shall calculate two- and three-body matrix 
elements between plane-wave external states. On the other hand, when calculating the 
matrix element of non-local operators, such as Lg, etc, we have to replace the simple 
plane waves by general wave packets. 

An important remark is in order. The total angular momentum operator in Eq. (IT^ . 
originates from the underlying 5*0(3) rotational symmetry of the NRQED Lagrangian. It 
has nothing to do with the perturbative expansion and must be exact. The gauge- invariant 
effective operators, as will be derived from the matching to the full theory, must agree with 
Eq. f lT2|) in the sum. This fact imposes the constraint S^*^ -|- + JJ^^ = Jnrqed, which 
in turn, leads to the sum rules among the above Wilson coefficients, as given explicitly in 
Ref. Q. 
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III. TWO-BODY MATCHING THROUGH SINGLE-ELECTRON STATES 



In this section, we consider matching the relativistic and NR angular momentum op- 
erators through single-electron states. This step alone cannot determine all the Wilson 
coefficients in Eqs. (fT3|) -( IT5|) . However, the calculation is useful to illustrate the general 
procedure, and shows simple physical concepts behind the matching. The result in this 
section is complementary to the three-body matching in the next one. 

We calculate the forward matrix elements of the angular momentum operators: {ep\Oj\ep) i 
{Oj)p,p, with Oj representing one of the relativistic operators S^, Lg and and \ep) is a 
free electron state with momentum p. We will perform non-relativistic reduction of their 
matrix elements in QED, and match them with those in NRQED. 

The spin operator is local and its contribution can be easily calculated in QED. At 
tree level, we have 

(Sg)p,p = (17) 

where u{p) is the usual Dirac spinor. The calculation of (Lg) and (J^) is more involved. The 
explicit presence of space coordinate x in the operator obscures the definition of the forward 
matrix element. Before going into explicit calculations, we address a basic question: what 
is the angular momentum for a free electron? 

The idea of the total angular momentum of a plane wave is ill defined. To see this, we 
first consider an off- forward matrix element of the operator Lg. After some algebra, 

(Lg)p/,p = J rf^x(p'|^^(x X 7r)^|p) 

= j c^^xe-^^'-^P'-P^x X i{p',p) 

= (27rf[2V,,53(p-p')]xf(p',p), (18) 

where i{p',p) = u{p')pu{p) at tree level. The derivative on the delta function is singular if 
one takes p = p'. 

A more realistic question would be to calculate the angular momentum of a distribution, 
or a wave packet |$) = J ^^$(p')|p'). The result will not only depend on the intrinsic 
properties of the electron itself, but also on the shape of the wave packet. Such relevance 
arises because the angular momentum operator is a non-local operator. In the momentum 
space, x can be Fourier transformed into the derivative over the conjugate momentum Vj,/. 
After integration by parts, the contribution to the total angular momentum is therefore 
separated into the two parts: the derivative on the matrix element i{p',p) and on the wave 
packet $(p'), 

(Lg)$ = [ rf3p$*(p)$(p) Uvp, X f(p',p)j d'pVp [$*(p)$(p)l X f(p,p).(19) 

J L Jp'=p Z J L J 

In the first term, [zVp' x i{p' ,p)]p'=p is the orbital angular momentum density intrinsic to 
the electron and contributes in the same manner as (Sg). The second term is precisely the 
angular momentum carried by the wave packet. Similar terms exist for the matrix elements 
of J^, as well as any non-local effective operator. In the following, we will denote the class 
of contribution like the first term in Eq. f|T9|) as the local contribution, and the second term 
as the non-local contribution for a matrix element between general wave packets. For local 
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operators like Sg, only the local contribution is present. For non-local operators, such as 
Lg and J-y etc, both local and non-local parts of the matrix elements from the full theory 
and the effective theory should match, independently. This will provide us with enough 
information to obtain the relevant Wilson coefficients. 

For non-relativistic reduction, we use the convention for a four- component spinor 



per 



2p0 



Uh 



(20) 



where Uh = \/2p^Q^ or 2p^ and p^ = \/p'^ + rn^- So we can rewrite the full theory 
matrix elements in terms of Uh and compare them with the effective theory. At tree level, 
the non-relativistic reduction of the QED matrix elements are 



\ '-''?/ local 
/T \(0) 



[(cr X p) X p] C(m" 



Aw? 



1 \u 



X p) X p] 0[m ^) 



f(p,p) 



(0) 

Lqjnon— local 



u\vuh, 



(21) 



This simple result can be used to examine the FW-tranformed and L^. By calculating 
the matrix elements of operators in Eq. ^ between two non-relativistic spinor wave packets, 
we find they indeed agree with Eq. fl2T]) . Therefore, we conclude 



(0) 
R 



(22) 



at the leading order. 






(b) 



(c) 



FIG. 1. One-loop corrections to (e|J|e) in QED. Here ® can be either Lg,Sg, or J-^: (a) for S^; 
(a)(b) for Lg; (c) for J^. Wave function renormalization diagrams, mass counterterms and the 
mirror diagrams are not shown explicitly. 



The one- loop Feynman diagrams for the electron two-body matrix element of Sg, Lg and 
Jy are shown in Fig. [H We use the dimensional regularization (DR) for both infrared (IR) 
and ultraviolet (UV) divergencies. The local contributions in the matrix elements are 

(S,)loL = ^"t{f + J^(-p)xp}., 
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where all IR divergences cancel. The non-local contributions are 



fb.P)S,non-local = (P. A, non-local = ^ (~3^ ~ ^ ~ y) ^^''^'^^ ' ^^^^ 

As a cross check, we sum up the one-loop results 

(S, + L, + J,)il, = 

f (P' P)L,^,non-local + f (P' P) J^^non-local = • (25) 

The vanishing sum of the first-order contributions is expected since the total angular mo- 
mentum is a conserved quantity. The algebra it respects dictates that the expectation value 
must be integer or half-integer and should not depend on expansion parameter a^^^. 

For the NRQED diagrams, we use dimensional regularization as well. Following the 
argument in Ref . [3] , in DR all loop diagrams in NRQED are zero due to scaleless integrals. 
The lack of scale is a natural consequence of the multi-pole expansion j20|. In the one- loop 
diagram containing ultrasoft photons in NRQED, the mass m decoupled and the relevant 
scales are IR cutoff and UV cutoffs in the effective theory. However, in DR both divergences 
are regulated by the extra dimension 4 — D and the only scale is /j,. Hence, all the NRQED 
loop integrals take the form 

. ((9,eff)(^)(/.) = (—--). (26) 

iTT \euv em J 

where wf''' is the zero-th order Wilson coefficient and (0i.cff)*'^'' is the first-order matrix 
element of effective operators, which is scale-independent. 
For a general matrix element in QED, we have 

{Oqeb^^ = ^(— + — + {A + B)\n^ + c] . (27) 
2n Veuv em 



where UV divergence reflects the renormalization property of the operator in the full theory. 
By imposing the matching condition, 

(C»nrqed)(')(/^) = {w,0,,,sf\fi) = ■ {0,,,Kf\fi) + ti;f ^ • (a,eff)^°H/^) 

- (—-—)+ wl'\f,){0,^sr- (28) 



27r Veuv em 

Since the full and effective theories have to reproduce the same IR divergences, matching 
yields, 

A,s = -B, ^«(;,)(O,,ff){0) = ^(A + i?)ln4 + C- (29) 

Threrefore, the Wilson coefficients can be obtained by subtracting all the 1/e poles while 
keeping logarithms and finite terms in full theory calculation. Effectively this is equivalent 
to "pulling up" the IR divergences to the UV divergences in the Wilson coefficients. Of 
course, this does not make sense since the Wilson coefficients should not depend on physics 
in the infrared. All the IR sensitivities that appear in the full theory calculation should be 
reproduced by the matrix elements of the effective operators. A serious problem resulting 
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from the above procedure is the obscure physical meaning of the cutoff scale yU in the di- 
mensional regularization. It yields the wrong scaling behavior of the matrix elements 

where Oj = Sq, hq and J^. To alleviate this confusion in QED, one option is to regularize 
the IR divergence using a different regulator, e.g., the photon mass. In the effective theory, 
the scaling behavior of the matrix elements also fails to satisfy the correct evolution equa- 
tionsfrom the naive matching in dimensional regularization. In Subsection III.F, we will 
render this problem by regularizing the UV divergence using a three-momenta cutoff A in 
NRQED, which is distinguished from the QED scale fi. We will also use the photon mass 
to regularize the IR divergences. 

At the moment, for the sake of simplicity, we use DR. 

In NRQED, the contributions to the electron two-body matrix elements can be readily 
obtained from the operators in Eqs. f|T3|) - f|T5|) . The /oca/ parts of matrix elements are 



2 4m2 



a X p) X p^Uh 



Uh 



(Lf )local = 4 + ^(O- X p) X p| 

(Jf )iocai = u[ j/^l + J^ia X p) X p| , (30) 
and the non-local parts of matrix elements are 

f(P,P)Lf, non-local = dRu[pUh , 
7°*^, non— local 

fRuipuh . (31) 

Comparing them with Eqs. (12T]) . (123|) and (124|) . we immediately arrive at some of the Wilson 
coefficients in Eq. ([I3D-([T5D in NRQED 



1 _|_ '^EM 1 I '-^EM 



27r ' " 27r ' 

" 27r V 3 m2 9 2n \ 3 9 Gtt ' 



•'^ 271 \3 m? 9 ) ^ 27r V3 m2 9 y Stt ^ ' 



These results agree with Eq. (7) of Ref. [Ml- Again, this result assumes dimensional regu- 



larization for IR and UV divergences in the effective theory calculations. 



IV. THREE-BODY MATCHING THROUGH ELECTRON SCATTERING IN 
BACKGROUND FIELD 

The two-body matching is insufficient to derive all the Wilson coefficients for the effective 
operators. For example, in two-body calculations, these operators involving at least one 
photon field give null result. To obtain the complete expansion, we have to consider more 
complicated matrix elements. In this section, we consider the processes with an electron 
interacting with an external electromagnetic field, through the angular momentum operators. 
Together with the results from the previous section, we will be able to determine the Wilson 
coefficients for all effective operators bilinear in the electron fields. 
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A. Tree Level 

We calculate the QED amplitude {epi\0 j\epKq) with an external vector potential A, as 
well as {epi\0 j\epA^^ with an external scalar potential . The in- and out-electron states 
are always taken to be on shell. We further assume electrons are non-relativistic and the 
virtual photons carry a small momentum q compared with the electron mass, i.e., <C m^. 




(a) (b) 

FIG. 2. QED tree diagrams for electron scattering in a photon background, coupled to angular 
momentum operators. 



In QED, the tree-level diagrams for the three-body matrix elements are shown in Fig. [2l 
Fig. [2] (a) is the direct contribution is from the x x A part in Lg. Fig. [2] (b) is an indirect 
contribution. We calculate the elastic scattering matrix element for on-shell electron = 
p>o _ p2_|_^2^ _ |p/|^ rpj_^g incident photon momentum is q'^ = (0, q), in which q = p' — p. 

Clearly the diagrams such as Fig. I2^b) will be divergent if we take p + q = p' since the 
intermediate state carries a momentum on mass shell. This is a pole term. During the 
matching, part of this pole term will be exactly reproduced by the non-relativistic electron 
propagator in NRQED. There is also a finite contribution though, due to the propagation of 
positron mode. According to Eq. (120|) . the virtual position propagation will be far off-shell 
and suppressed by 0{m~^). Effectively, the propagator of Fig. EJ^b) will shrink to a point 
and result in a contribution from a local effective operators like those in Fig. |2]^a). The 
physics beyond the scale m shall be be in the corresponding Wilson coefficients. 

To regulate the above pole term, we introduce a small off-shellness to the intermediate 
electron by assigning the photon momentum q'^ a time-like component, q'^ — )■ q). Mean- 
while we rewrite the intermediate angular momentum operators as (9 — ?■ Oe^''°^° to maintain 
the four-momentum conservation. Diagrammatically, we have an additional small momen- 
tum [q^, 0) flowing in from the photon and out from the angular momentum operators. 
Eventually, we will take g° — ?■ 0. 

For the matrix elements, we separate them into the 1/q^ pole term and the regular part 
in which represent the contribution from non-local electron propagation and the local 
interactions in NRQED, respectively. In practice, it is enough to match only the local part 
of the matrix elements between QED and NRQED. Matching the terms carrying poles 
can be used as a consistency check of our results. 

In the presence of g", the denominator of the intermediate electron propagator in Fig. MJo) 



can be Taylor expanded as [11 



1 1 1 g° p2 



(p + g)2 — 2mq^ Am? 9>m? Am^q^ 

where we have used the on-shell kinematics for the electron p"^ = p'"^ = m^, p' = p + q and 
p'O = pO ^ 777, -|- p2/2m and make subsequent non-relativistic expansion. Later on we shall 
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make the same expansion in powers of g° in the next-to-leading-order calculations as well. 

For non-local operators, we choose x — t- zVg, i.e. the derivative with respect to the 
external momentum of the photon. This has several advantages. First, the result will 
be symmetric with respect to the incoming and outgoing electron momenta. Second, we 
need not consider the derivative on the photon polarization which does not depend on q. 
This greatly simplifies the calculation. Finally, the three-momentum conservation relation 
p' = p + q should be understood after the replacement x — )■ zVg has been made. 

The tree-level non-pole contributions from Fig. [2] in QED are 



(0) 
local 



/T \(0) 

\-^q/ local 



-^ui [iq X A + (P ■ A)cT - (P ■ a) A] uh , 
[-iq X A - (P ■ A)cr + (P ■ cr)A] Uh 



i{p',p,p -p) 



(0) 

Lq, non-local 



1 - 



|q| 



5m^ 



+ 



5m^ 



(qxP: 



'-eA)uh 



(34) 



in which P = p' + p and q = p' — p. Terms of order 0{m ^) have been suppressed. 



p+q 



p+q 



(a) (b) (c) 

FIG. 3. corresponding NRQED tree diagrams as in Fig. 2 



Next, we consider the tree diagrams in NRQED as shown in Fig. |3l The electron propa- 
gator in Fig. |3]^b),(c) 



Ep — 



(p+q)^ 
2m 



(35) 



is also divergent for on-shell electron external states. By introducing a small non-zero 



the fermion propagator will be regularized to 



Since there is no angular momentum 



effective operator containing time derivatives, the pole can only be canceled by Dirac 
and spin-orbit terms in the NRQED Lagrangian involving the electric field strength E (since 
E* = —F^^ ~ 2(5° A* — gM'') contains positive powers of g°), resulting in local matrix 
elements. 



As a concrete example, we calculate the local part of the matrix element (e| J ci^x^/^^x 



X 
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TzipleA)^^^'^'^^ as shown in Fig. [3|, 



(e|x X 7r|eA) 



(0) 
local 



iVg X ul < i-eA) 



+ |^g°([(2p'-q).Axa]- 



+ ie 



+ 



-q^ + ie 



5m^ 



q ■ (qA° - g°A) 



ip 



le 



[(2p + q).Axo-] 



zp' 



p -p 



lecD 



iecf) 
8m^ 



g°A 



P + P 



ec5 

r 

i2 



g°(A X 0-) 



Uh 



q X A 



ecs 



5m^ 



P X fcr X A) 



(36) 



In the above calculation, we have included both cases that the NRQED interaction happens 
before and after the angular momentum operator, along the momentum flow. As promised 
above, we only need to keep the non-pole terms [regular as g° — ?■ 0] and suppress all the 
poles. 

Following the same procedure, we have calculate non-pole contributions from the local 
and non-local matrix elements for all the effective operators in Eqs. f lT5]) - ffT^ . The results 
are listed in Table I. 



operator 


local 


non-local 




-^csuliP X A) X auh 






-8^ul [cD(«q X A) 
+csP X (cr X A)] Uh 


8^4 [Cz5(q- A)q 
+icsCT ■ (P X A)q] Uh 


i/j^ [x X (o- X E)] V 


2uliA^a)uh 


-iul{A°cr X q)^;, 


Ip'^ [{a X Tt) X TT 

—TV X (cr X tt)] ip 


eu[ [P X (cr X A) 
+A X (cr X P)] Uh 






m^q X Auh 




X X [B X Viip'^ip)] 


mj^q X Auh 


ul [q2 A - (q • A)q] Uh 


X X ^/jt [B X (cr X ^)] -0 


njjA X (cr X P)uh 


—iu\ [cr • (q X P)A 
+cr • (P X A)q] Uh 


X X [cr X (B X ^)] ip 


ul [(a • P)A + (A • P)a] Uh 


iu[ [(q • P)cr x A 
-(A • P)cr X q\uh 



TABLE I. Local and non-local contributions to (elO |eB)Trcei only local terms are shown 



From Eqs. (fT3|) . ( !T^ and (IMj) . we get the leading order Wilson coefficients up to 0{1) 

= -df^ = -df = rfg) = df = 1 . (37) 

They agree with the non-relativistic reductions through the FW transformation in Eq. (Q. 
The other Wilson coefficients in S^*^ and L^''^ are of order (9(ag^) or higher. They will show 
up in one-loop matching as described in the next subsection. 
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B. One-Loop Matching 



In this subsection, we calculate the one-loop radiative corrections to angular momentum 
operators in QED and NRQED. To find the complete set of Feynman diagrams, we first 
draw the one-loop QED diagrams and insert angular momentum operators in all possible 
ways. The corresponding Feynman diagrams in QED are shown in Fig. HI Here we have no 
diagram with on the external photon line because this yields only pole contribution. The 
photon pole contribution is related to matching in the photon state, which will be done in 
the next subsection. 




(i) S) (k) 



FIG. 4. One-loop contributions to (e|J|e7) in QED. Here (S" can be either L|j,Sg, or J^: (a)-(i) for 
Lg; (b)(d)(g)(i) for S^; (j)(k) for J^. Wave function renormalization diagrams, mass counter terms 
and the mirror diagrams are not shown explicitly. 



The loop diagrams are considerably more complicated than the tree-level case. This is 
because all loops can be Taylor expanded in g° and contribute to the local matrix elements 
in a similar manner as Eq. f l36p . Again only the regular part of matrix elements in the limit 
g'^ — )■ is taken into account. 

After lengthy calculations we find, for the vector external potential A case, the local part 
of matrix element reads 

(S,)Lii = ^ A«I [7^q X A + (P . A)^ + . A)P - 2(P ■ ^) A] u, , (38) 



(J 

V-^g/ local 



27r 



8 In 



em 



50 

y 



zq X A 



fP ■ A)cr - 2fcr ■ A)P 



fP ■ (t)A 



Uh , 
(39) 



(1) 



7 /local 



2tt Ath? 



-— -81n^-^) zqx A-^(P-A)^ + (a-A)P-|(P-a)A 

eiR m"^ 3/3 3 



Uh , 

(40) 
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and the non-local part 



f(p',p,p'-p)2,non-Iocal 



3euv 



(1) 

J^, noil— local 



"^EM ^ t 

27r 4m2 

--z[cr- (P 

3 ^ ^ 



--(q- A)q + 



"^EM ^ t 

27r 4m2 



+ -2[o--(Px A)]q- 



4 17\ 2 

+ - In + — ) (4m2)A - -^(P ■ A)(cr x q) 





3euv 3 


2 


4 




3eiR 


4 




3euv 


3 m-^ 



9 



In — 



31 
18 



z[cr • (q X P)]A 



2 In At - - 

6 



17 



|q| A 

2 

L + -i 
3 



(41) 



n At + — W o- ■ q X P A 

3euv 3 m2 IP ' L Ji 



+ g(q- A)q- 



|q| A 



3euv 3 6 
For the scalar potential case, the local part of matrix elements are 



Uh ■ 



(42) 



(1) ^ Q^EM e t 
^""'^^ 27r 4m2 ^ 



4m 



(1) 

local 



27r 4m2 



A^cr 

3 



and the non-local part 



i{p',p,p' -p)l , non-local 



27r 4m2 



Uh + 0{m 



2"^ ,0-/ X 
-— A°z(cr X q) 



(43) 



i{p',P,p' -p) 



_ Q'em ^ t 
J^, non-local 27]- 



(1) 



2m 



y4°i(cr X q) 



Uh 
Uh ■ 



(44) 



Again we sum up all local and non-local matrix elements of the angular momentum and 



get 

(e|Sg 
(e|S„ 



Lg + J^leA)!^^^! 

Lg + J7|eA)|j2i-local 



e\Sg + Lg + 3,\eA')[ll, 



3^\eA^) 



0\(1) 

non-local 



27r 4m2' 
. 



0, 

2^[^-(qxP)]A-^f— + lni^) |qpA 
3 VejR 



The reason for the null results is again due to the conservation of the total angular mo- 
mentum. The non-zero first-order result for {e\3-^-EiQED\e-^)^orhitai non- vanishing 
one-loop contributions of the electromagnetic form factors Fi{q^) and F2{q^). Once we in- 
clude the diagrams with the operator placed on the external photon line, the total orbital 
contribution to the matrix element (e| J7|eA)^],|jj^jj[ would also yield zero. 

For the matrix elements in NRQED, we only need to calculate the tree-level diagrams for 
the effective operators in dimensional regularization. By applying the matching conditions 



16 



according to Eq. (fT6|) together with the results in Table I, all the Wilson coefficients for 
effective operators bilinear in quark fields can be determined in the MS scheme, 

1 I ^EM „ -r I "em „ 1 I ' "em 

J- + -T, — ! Ob = i + 



^ 27r V 3 m2 9 / ' 27r V 3 9 / ' Gtt 

ciz. = l + ^f-41n4--V ci5 = l + ^f--ln4 + - 
27r V m2 Sy ' 27r V 3 m2 9 



EM 



. "em /4 /i^ i7\ Q,^^^ /4 11\ 

^^ = ^U^^J' ^"^^W^'^^J' ^^^^ 
r "em /4 /i^ 1\ / 4 u2 13 





) fcr 


"em 

27r 




fs = - 


"em 

27r 


"em 

37r ' 


fB = 


"em I 

27r 



As it can be easily verified, the coefficients satisfy the sum rules in Eq. (20). 

These results reproduce Eq. (8) of Ref. except for coefficients dn, ds and fs- We 
will properly separate the UV and IR divergences in these three coefficients when discussing 
the renormalization group evolution in Sec. VI. 



V. PHOTON SECTOR 

In this section, we calculate the forward matrix elements (71017) up to the order 0{aj^^^) 
with off-shell photons carrying a momentum q. The tree-level matching is straightforward 
and we have 

/f = 1, 4°) = 0. (47) 





*W\A/ 




FIG. 5. One-loop contributions to (7|Jj7) in QED. Mirror diagrams are not shown, (a) for Sq; 
(a)(b) for Lq] (c) for J^. 



The relevant QED Feynman diagrams at one-loop are shown in Fig. O While the op- 
erators Sq and Lq can be inserted to the fermion loop of the vacuum polarization diagram 
as in Fig. [5]^a) and (b), the one-loop contribution from only arise in Fig. ^c). We need 
no new effective operator to match the matrix element of J^, simply because the vacuum 
polarization is already reflected by the di- and (i2-terms in the effective Lagrangian. Thus 

= , (48) 
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and there are no further corrections to at 0(a*^ m ^) order as aheady indicated in 
Eq. ([ISD. 

The calculation of the one-loop matrix elements for Sg and Lg in QED yields 



(7|S,|7)!oLi 

(7 1 1 7) local 



(1) _ "em 



+ 



[(A X A*)q\' - (A X q)g2A*° - (q x A*)gM°] , (49) 
2-71 \rn? J 3euv 



^ TT^ [(A X A*)g° + (A X q)A*° + 2(q X A*)A°] 



(1) 

Lo, non— local 



, ^ [(A X A*)g%2 + 3(A x q)^^^*" + 13(q x A*)q^A'^ 
In 15 

+ 4(A X q)(g ■ + 4(q x A*)(g ■ A)q^] , 

2 \ ^ 2 

[-2(A ■ A*)q\ + (g ■ A)q^A* + (g ■ A*)g°A] 



(50) 



a 



+ 



/r 

27r \rn? ) 3euv 
1 2 



[-4(yl ■ yl*)g°g2q + . A)q^q^A* + (g ■ A*)q^q^A 



271 m? 15 

+ (g ■ A)A*\^ci + (g ■ A*)AW - {q^iA^^A* + A*°A) 
+ 2(g-yl)(g-A*)g°q 



(51) 



where A and A* represent the incoming and outgoing photon fields, respectively. 



operator 


local 


non- /oca/ 


X X (E X B) 


i [(A X A*)gO + (A X q)^*o 
+2(q X A*)A0] 


-2{A ■ A*)gOq 
+(g- A)g°A* + (g-^*)gOA 


{d^FP^)Fp^ 


-2z [gOg^A X A* 
-{q ■ A*)q^A X q - (g • A)gOq x A*] 




{d^Ff^)Fp^ 


-i [(g2^0 - q ■ V)q X A* 
+ {q^A*^ - q ■ A*q°)A x q] 




(UF^P)Fp^x' 


-i [4g2^0q X A* - q^A*'^A x q 
+q^q'^A X A* -2q- Ag°q x A*] 


-g2(g. A*)(gOA + AOq) 


{d^FP''){d^Fp^)x^ 


4zg°(g • A)q X A* 





TABLE II. Local and non-local contributions to off-shell matrix element (7|0'^'^|7)Tree- 



In NRQED, since the energy fluctuation cannot excite the virtual electron-position pair, 
there are no loop corrections to the matrix elements of the effective operators. All the 
loops in QED will be encoded in the corresponding Wilson coefficients. All the tree-level 
matrix elements of the effective operators in Eqs. f|T3|) . f|T^ and f|T5|) are summarized Table 
II. Matching the matrix elements between QED and NRQED, it is easy to find the follow 
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combinations of effective operators 

2 



ij3 



15 



15 



15 



(52) 



A little algebra yields the operator relations 

{d°FP')Fj -^ExE-BxB, 
{d'F^'^)F^ (V ■ E)B - (B ■ V)E 
UF^PFjx' -X X (DE X B) , 



pry 



X V)E'' - B"(x X V)B'' 



(53) 



Together with Eqs. (fT3l) and (fT4|) . we extract the Wilson coefficients in the photon sector in 
the MS scheme 



a. 



a. 



^71 



127r 



■^72 



67r 



/7 



^In— 

111 , 

1 . 



eOvr 



^72 



^73 



157r 



"74 



157r 



(54) 



VI. DISTINGUISHING QED AND NRQED UV CUTOFFS 

In the previous sections, we have obtained the non-relativistic reduction of the angular 
momentum operators with the main results in Eqs. (|T3|) -( fT5l) and Eqs. ( 146|) . ( l54l) . We have 
used the dimensional regulation for both UV and IR divergences in QED and NRQED, and 
the Wilson coefficients depend on a single energy scale yU.. However, there is no need that the 
NRQED must use the same UV regularization. In fact, in many bound state calculations, it 
is better to use three-momentum cut-off for the UV divergences and the photon mass for IR 
divergences. In this section, we will try to redo the calculations in this new regularization 
scheme. Before that, let's consider a problem that arises from DR for all divergences. 

We expect the UV behaviors of QED operators are completely included in the Wilson 
coefficients in NRQED. In another word, the matrix elements in NRQED must satisfy the 
same renormalization group (RG) evolution equations 

dt \3f) 2vr V I -\ + \) VJf ; ' 

where t = In jj^ /rn? and J^*^ = S^^ + L^^ represents the total angular momentum carried by 
the electron. The additional | in the evolution equation stems from the scale dependence 
of the redefined photon fields in the effective theory, the di term. 

Our effective operators in Eqs. f lT3|l - f|T5|) with Wilson coefficients listed in Eqs. fH6|) and 
f l5^ . actually, fail to satisfy the desired evolutions. The reason for this, as already discussed 
in subsection III, is due to DR. In the QED calculations, we have not distinguish the 



(55) 
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dependence from UV or IR divergences. The IR divergences should be captured by the 
effective theory and has nothing to do with the RG flow. On the other hand, when matching 
to NRQED, we also have not distinguished the UV cutoff dependence from that in QED. 
By introducing the three-momentum cutoff A in NRQED calculations, we can restore the 
correct RG evolutions. In addition. We let the photon have a small mass A in both QED and 
NRQED to regulate the infrared. According to Ref. [13] the Wilson coefficients in NRQED 
Lagrangian with three- momentum cutoff A as UV regulator are 

c. = c, = l, c, = l + ^, cz. = l-^An^-?V cs-^ ' 



d, = l + ^\n^, d, = ^. (56) 

The separation of two scales fi and A is obvious. 

The full theory amplitude with the photon mass regulator can be translated by making 
the replacement in Eqs. f l38|) - fH2l) 

— + ^41 ->21nA . (57) 

In the effective theory, the loop diagrams are no longer vanishing if we use A and A to regulate 
UV and IR divergences. We choose to work in the Coulomb gauge to do the calculation since 
the full result is gauge- invariant. All the diagrams under consideration in NRQED are listed 
in Fig. [6l The non-relativistic electron propagator picks up only the pole at = \/p'^ + rn^- 

^^^^""^^^^ E, + k^ _^M^ + ,, = FT^+2m(p|ze)2+---- 

This can be understood as a multipole expansion in inverse powers of m. 
The one-loop effective theory matrix elements now read: 

{\ )iocai - ^ V A 6 J ^ ' 

/Jefr^(l) _^l!^^(o)An^_5\ t„xAM. 

(Lf )itL,ocal = T - " ^ "^"^"'^"'^ • (^'^ 

The one-loop amplitude (S^^)iocai and (J^j,''^)non-iocai as well as the diagrams involving the 
scalar potential vanish at higher order 0{m~'^) in NRQED. 
Now the matching condition Eq. ( !T6|) should be rewritten as 

{J{fi))QEB = {r\fl,A))^KQED. (60) 

The sum rules of the Wilson coefficients in Eq. (??) do not change. However, all the Wilson 
coefficients in the effective Lagrangian are defined in the presence of new regulators, similar 
to Eq. dM]). 



20 




(i) (i) (k) 



FIG. 6. One-loop corrections to (e|L^f^|eA), (a)-(f), and {e\3f\eA), (g)-(k), in NRQED. In the di- 
agram the fihed box represents the vertex interaction from term — ^ in the NRQED Lagrangian. 
In the diagram the fiUed box represents the vertex interaction from term — ^ in the NRQED 
Lagrangian and the dot represents the vertex. The dashed hne is for the Coulomb photon 
propagator while the wavy line is for the transverse photon propagator. Wavefunction renormal- 
ization diagrams, mass counterterms and the mirror diagrams are not shown explicitly. 



With the new matching conditions, we find yet another set of Wilson coefficients depend- 
ing on both three-momentum cutoff A in NRQED and the QED cutoff fi. We only list the 
ones different from those in Eq. P6|) . and denote them with an asterisk. 



D 



rl* 



2tx 



a. 



ih4 

3 m-^ 



+ 



2tx 



-16 In 



m 
2A 



16 m 
— In — - 
3 2A 



23 



37 



B 



2ti 



16 In 



m 
2A 



86 



(61) 



With the new coefficients, both J^^ and J^ ^ now satisfy the RG equations, Eq. (l55l) . They 
now agree with those in Eq. (8) in Ref. 14 . 



VII. CONCLUSIONS 

In this paper, we have established systematically the relativistic angular momentum com- 
ponents in the framework of NRQED up to the order a^^/m^. Such effective operators can 
be used in the computation of higher-order contributions to spin/orbital angular momentum 
in non-relativistic bound states, such as the hydrogen atom. Further extensions to the ef- 
fective operators in NRQCD can be readily obtained following the same procedure, and will 
be useful for studying the strong-interaction bound states such as the heavy quarkonium 
systems. 
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Appendix A: Generalized Noether Currents with Higher Derivatives 



If the effective theory Lagrangian contains higher derivatives terms the conventional equa- 
tion of motion and Noether currents, derived from L = C{(j), d'^(j)), have to be modified. In 
this appendix, following the procedures in Ref. 23|, we have generalized the formulas to a 
Lagrangian with second-order derivative C = C{(j),d^(j),d'^d'^(f)), which will be relevant to 
deriving the angular momentum in NRQED, Eq. f[T^ . We list the results here. 

We use the following notations: 



dC 



(Al) 



The generalized equation of motion is 

^ . (A2) 

The generalized, Belinfante-improved, symmetric energy-momentum tensor T^'^ has the form 

- dj^''" , (A3) 



where f'^'^'^ is defined as 
1 
2 



/ 



fiau 



+ ^[Co.rd,^p){S"^r - (/i o a) + (/i ^ z/) 



(A4) 



where (5'"^)'^'' is the generator of the Lorentz group. Finally, the generalized angular mo- 
mentum tensor can be obtained through 



(A5) 



As an application to NRQED, for the higher derivative operator ^F^y,\3F^^^ in the La- 
grangian, we obtain the corresponding angular momentum operators 

2d2 



Jj3 



2UF^PFjx' + d^FP'^&Fp^x' + 2{d^FP')F- 



-X X (DE X B) + E"(x X V)E" - B'^fx x V)B" + ExE-BxB 



(A6) 



which corresponds to the terms proportional to d2 in Eq. (fT2|) . 

The other terms in £nrqed contains at most one derivative, so the remaining parts of 
angular momentum effective operators can be readily obtained. 
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